Dynamical instability, Chaos, and Bloch oscillations of Bose-Einstein condensates in 

tilted optical lattices 



O 
O 

(N 

o 

Q 

00 



Andrey R. Kolovsky 
Max- Planck- Institut fiir Physik komplexer Systeme, D-01187 Dresden and 
Kirensky Institute of Physics, Ru-660036 Krasnoyarsk. 
(Dated: February 2, 2008) 

We study the Bloch dynamics of a Bose-Einstein condensate of cold atoms by using the formalism 
of the discrete nonlinear Schrodinger equation. Depending on the static force magnitudes the system 
is shown to exhibit two qualitatively different regimes of Bloch oscillations - exponential decay for 
the static force magnitude less than some critical value (defined by the condensate density) and 
quasiperiodic oscillations in the opposite case. The relation of these regimes to the onset of chaos 
in the system is discussed. 
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I. INTRODUCTION 

Recently much attention has been payed to Bloch oscil- 
lations (BO) of Bose-Einstein condensate (BEG) of cold 
atoms in the optical lattices [illlllillli[llalili3- 

It was found that, unlike the case of a dilute atomic 
gas, where the atom-atom interactions can be neglected, 
BO of BEG rapidly decay. A plausible explanation of 
this phenomenon is suggested by the dynamical insta- 
bility of the Gross-Pitaevskii equation in the presence 
of periodic potential (see papers 0, Q and references 
therein). In this paper we revisit this approach, assuming 
for simplicity the relatively deep optical lattices, where 
the Gross-Pitaevskii equation can be approximated by 
the discrete nonlinear Schrodinger equation (DNLSE). It 
is shown below that a strong static force can actually 
suppress the dynamical instability of the DNLSE and, as 
a consequence, BOs become not-decaying. We indicate 
a condition on the critical magnitude of the static force 
separating the decaying and not-decaying regimes of BO 
which, for a fixed static force, can be also formulated as 
a condition on the critical density of BEG. This result is 
well supported by the numerical simulation of the system 
dynamics, performed for the realistic parameters of the 
present-days laboratory experiments with cold atoms. 



II. DYNAMICAL INSTABILITY 

According to DNLSE the (mean field) Bloch dynamics 
of a BEG is described by the system of coupled nonlinear 
equations 



ihai 



{ai+i+ai^i)+g\ai\^ai-^dFlai . (1) 



In Eq. ai{t) is the complex amplitude of a BEG of 
atoms localised in the Zth well of the optical potential, J 
is the hopping or tunnelling matrix elements (who's value 
is uniquely defined by the lattice depth) , d the lattice pe- 
riod, F magnitude of the static force, and g the nonlinear 
parameter given by the product of the microscopic inter- 
action constant W and the filling factor n (mean number 



of atoms per lattice site), g — nW. For a nonuniform 
density n will mean the peak density. 

Let us first consider the case of uniform initial condi- 
tions a/(0) = 1. (Note that we do not normalise J^i 
to unity.) Subsequently using the gauge and Fourier 
transformation, b^, = L^^/^ ^^^-^ exp(iK^ — iijjBlt)ai, 
Eq. 1^ can be written in the form, 



ihbK. = — J cos{dK — LL!Bt)b^ 



(2) 
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where k is the quasimomentum {—n/d < n < Tr/d), and 
lub = dF/h Bloch frequency. For the given initial condi- 
tion Eq. has trivial solution 

bo{t) = cxp (^i-j^ sm{ujBt) - i^t^ y/Z , bj^^oit) = , 

(3) 

which is no other than the celebrated BO, where the ki- 
netic and Stark energies of the system oscillate according 
to the cosine law with the Bloch period Tg = 2tt/ojb = 
h/dF . However, for g the solution (PJ can be unsta- 
ble with respect to small perturbations. Following the 
standard approach we linearise Eq. around the solu- 
tion (PI, which leads to the system of linear equations. 



ihb+K = — J cos{dK — uiBt)b+K 
ihb^fi = — J cos{dK + ujB't)b-K 
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-blb\^ ,(4) 



where &±k(0) are arbitrary small. Substituting here 60 (i) 
from Eq. Q and integrating Eq. in time, we have 
b±n{t) ~ eyi'p{vt)b±^{Qi), with v given by the logarithm 
of the maximal eigenvalue of the Floquet matrix 



U — exp 



J{t) = exp 



Tb 



1 f{t) 



dt 



(5) 



2J 

i— — [1 — cos(dK)] sin(wBi) 
dF 
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FIG. 1: Increment of the dynamical (modulation) instability 
V as the function of the static force magnitude F and the 
quasimomentum k for g — 0.4. The lattice period d and 
the tunnelling constant J are set to unity. The inset shows 
u — v[F) for K = n /d (edge of the Brillouin zone). 



where the hat over the exponent denotes the time order- 
ing. It is easy to see that the matrix U is parametrised 
by the product dn, the ratio g/dF, and the ratio J/dF. 
Thus, without lost of generaUty, we can set two parame- 
ters (J and d in what follows) to unity. A typical depen- 
dence of the increment v on the quasimomentum n and 
the static force magnitude F is depicted in Fig. It is 
seen in the figure that there is a critical value Fcr ~ 1.2 
above which v = Q and, hence, the solution Q is sta- 
ble. Scanning over different 5, we found Fcr ~ 3.g for 
0.1<c/<l|l. 

Suppression of the dynamical instability is easy to 
check numerically. These numerical simulations, per- 
formed for different lattice size L and periodic boundary 
conditions, undoubtly indicate existence of the critical F. 
Moreover, this is also valid for the non- uniform initial 
conditions jajp — exp(— /^/cr'^), realized in the experi- 
ment. Although in this case and for F > Fcr BO slowly 
decay (with subsequent revivals, see Fig. 0| below), this 
decay is due to the trivial dephasing [g and fundamen- 
tally differs from the exponentially fast destruction of BO 
at < Fcr, caused by the dynamical instability. 



III. TWO-SITE MODEL 

One obtains a useful insight in physics of the discussed 
phenomenon by considering the limiting case of the lat- 
tice with only two sites, 

2 

H{t) - -Jcos{ujBt){a;ai + 0^03) + f I] hi" - (6) 



FIG. 2: Phase portrait (stroboscopic map) of the system ^ 
for g — 0.4 and F = 2.0 (upper left panel), F = 1.3 (upper 
right), F — 0.7 (lower left), and F = 0.5 (lower right panel). 

where |aip4-|a2p = 2 is the integral of the motion. Using 
the substitution ai.2 = ^J1I\,i exp(i0i^2), I = h ^ ^2, 
9 = 9i — 62, one maps (0 to the periodically driven 
pendulum, 

H{t) = gl^ - Jy/l - P cos(wBt) cos 9 , (7) 

with the angular momentum / resticted to the interval 
l-^l < 1 11]- The stroboscopic map of the pendulum O 
is depicted in Fig.|5]for different values of the static force 
F = 2.0, 1.3, 0.7, 0.5. The periodic solution ®, which 
is of our particular interest, corresponds in Fig. |21to the 
stationary point {1,9) ~ (0,0). As F is decreased, the 
stability island surrounding this elliptic point shrinks to 
zero and at F w 0.9 (see inset in Fig. ^) it bifurcates 
into the hyperbolic point. We also note that for F < 
0.9 the phase space of the system is dominated by the 
chaotic component. In this sense the onset of dynamical 
instability at F « 0.9 corresponds to the condition of 
developed chaos [l^ . 

Next we discuss the relation between the onset of 
chaos and decoherence (depletion) of BEG. To address 
this problem we need to treat the system microscopi- 
cally, where for the discussed two-site model the micro- 
scopic Hamiltonian is given by the time-dependent Bose- 
Hubbard model, 

i W 

H{t) ^ - J cos{ujBt)ialdi+h.c.) + — 2_^hi{ni-l) . (8) 

1=1 

As an overall characteristic of decoherence we consider 
the trace of the squared density matrix, p^(t) ~ Tr[p^(i)], 
where 

pi,m{t) ^ N-^^it)\dlam\-^it)) (9) 

is the one-particle density matrix. (We recall that for a 
coherent evolution p^(i) = 1.) In Eq. Q {^(t)) is the to- 
tal A^-particle wave function of the system {N = 100—400 
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FIG. 3: Upper panel; Dynamics of p^[t) = Tr[p^(t)] for the 
system Q with n = 50 and g = nW — 0.4. Two curves 
corresponds to F = 1.3 [p'^{t) ~ 1] and F — 0.7 [decay of 
p^(t)]. The lower panel is the 'classical' simulation of the 
decoherence process. 
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FIG. 4; Bloch oscillations of the mean atomic momentum for 
f = 1.3 > Fcr (a) and F = 0.7 < Fcr (b). [Average over 
100 different realizations of the initial conditions IllH with 
n = 200.] The width of the atomic array is o" = 26.5, the 
nonlinear parameter g = nW = 0.4. 



in the numerical simulations) and as an initial condition 
|\l/(0)) we choose the ground state of © for F ~ 0, which 
is approximately given by the product of N Bloch waves 
with zero qusimomentum. It is found that when the cen- 
tral stability island in Fig. |21 is large enough in compar- 
ison with the quantum of the phase space volume 27r?i 
[T^ decoherence is negligible. On the contrary, when 
the periodic point (/, 9) = (0, 0) is hyperboHc, we ob- 
serve rapid decoherence of the system (see upper panel 
in Fig.|3Jl. Note that the rate of decoherence can be well 
estimated by using the classical (i.e., mean field) dynam- 
ics and considering an ensemble of initial conditions scat- 
tered around the periodic point. Then the one-partical 
density matrix @ is given by the correlation matrix, 

pi,mit) = L-\{aUt)aUm , (10) 

where the double-angle brackets denote the average over 
an ensemble of the initial conditions. The characteristic 
width of the distribution function for the action and angle 
variables is obviously given by the effective Planck's con- 
stant or, equivalently, by the quantum fluctuations of the 

number of atoms in one well, A/^ ~ h An^ ^ 1/n, and 
phase fluctuations, A6^ ^ Xj Av? ^ 1/n 15]. For n = 50 
the classical dynamics of [t) is depicted in the low panel 
of Fig. 13 By comparing both panels one concludes that 
decoherence of the quantum system is actually due 
to the chaotic dynamics of its classical counterpart JBJ. 

IV. BLOCH OSCILLATIONS 

Although it looks problematic to extend the above 
phase-space analysis of Sec. IIIII for L ^ 1, some con- 



jectures are still possible. Indeed, using the polar pre- 
sentation for complex amplitudes a/ , one maps the L-site 
system into a system of L driven coupled pendula and 
there is no doubt that this system would be generally 
chaotic . On the other hand, absence of the dynamical 
instability for F > Fcr indicates the presence of the stable 
periodic trajectory in the multi-dimensional phase space 
of the system corresponding to the stable BO. The crucial 
point is, however, whether the stability island surround- 
ing this periodic trajectory is large enough to support the 
quantum state. To answer this question we numerically 
solve the DNLSE for randomly varied initial conditions. 



91=^ (4n|az|2) , ai = \ai\exp{i9i) (11) 

where f is a random Gaussian variable with = 1 , which 
accounts for the quantum fluctuations. The results are 
presented in Fig. 01 which shows dynamics of the mean 
atomic momentum for g = 0.4, n = 200, a = 25.6, 
F = 1.3 (upper panel) and F = 0.7 (lower panel). It 
is seen that the stable regime of BO is also stable against 
the quantum fluctuations. This validates our conclusion 
about two qualitatively different regimes of BO. 

It interesting to estimate critical values of the param- 
eters for a typical laboratory experiment. Taking as an 
example the recent experiment Q with ^^Rb atoms in 
a vertically oriented quasi one-dimensional lattice of the 
depth 4 recoil energies, one has dF = 0.28, J = 0.17, and 
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W = 0.71 • 10~^ recoil energies, respectively. Then the 
condition on the critical magnitude of the static force is 
formulated as a condition on the critical density of the 
condensate and corresponds to ficr = 120 atoms per lat- 
tice site. Since the peak density in Ref. is larger than 
•her, BOs should decay, which agrees with the experimen- 
tal finding. 

V. CONCLUSION 

We have studied the dynamics induced by a static force 
of a BEC of cold atoms in an optical lattice. Depending 
on the static force magnitude F (or, equivalently, on the 
density of the condensate n) the system is shown to have 
two qualitatively different regimes of BO - exponential 
decay of oscillations for F < Fcr{n) and 'quasiperiodic' 
oscillations for F > Fcr{n). It is argued in the paper 
that the former regime reflects the chaotic dynamics of 
the system, where the transition to chaos coincides with 
onset of the dynamical instability in DNLSE. 

It is worth of noting that the above results are obtained 
within the mean field approach. Because this approach 
assumes the limit n oo, Wn = g/n — > 0, for any 
finite n the classical (macroscopic) dynamics of the sys- 
tem, governed by DNLSE, deviates from the quantum 
(microscopic) dynamics, governed by the Bose-Hubbard 
model. It is a challenge both for theory and experiment 
to establish the relation between the classical and quan- 



tum results, depending on the parameter h = 1/n. We 
would also like to stress that for considering the problem 
of quantum-to-classical correspondence one has to aver- 
age the solution of DNLSE over an appropriate ensemble 
of the initial conditions. Without this additional proce- 
dure, the solutions of the DNLSE may have nothing to 
do with the solutions of the Bose-Hubbard model. 

In this present paper we analyse the problem of 
quantum-to-classical correspondence for BO by using a 
two-site model. Both the similarities and discrepancies 
were found. In particular, we show that the classically 
chaotic dynamics of the system for F < Fcr is responsible 
for decoherence of the system, which is the deep reason 
for the decay of BOs. At the same time, the rate of deco- 
herence, obtained by the means of the DNLSE, appears 
to be correct only for the short-time dynamics. 

To conclude, we briefly discuss the relation of the above 
results with those of our recent papers H, 0, , devoted 
to the Bloch dynamics of the system in the deep quantum 
regime h ^ 1/n ~ 1 • Clearly, in this case one cannot 
appeal to the results of the classical (DNLSE) analysis 
and the system should be treated microscopically. Re- 
markably, that this microscopical analysis also reveals 
two regimes of BO - regular (not-decaying) oscillations 
for a strong forcing and rapid dec ay ( decoherence) of 
BO in the case of a weak static force [l3|. 

The author acknowledges discussions with J. Brand, 
A. Buchleitner, S. Flach, and S. Sinha. 
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